ROBUST AUXILIARY VECTOR FILTERING WITH CONSTRAINED CONSTANT MODULUS 

DESIGN FOR BEAMFORMING 



Lei Wang and Rodrigo C. de Lamare 

Communications Research Group, Department of Electronics, University of York, YOlO 5DD, UK. 

Email: {lw5 17,rcdl500} @ohm.york.ac.uk 



ABSTRACT 

This paper proposes an auxiliary vector filtering (AVF) algorithm 
based on a constrained constant modulus (CCM) design for robust 
adaptive beamforming. This scheme provides an efficient way to 
deal with filters with a large number of elements. The proposed 
beamfoiTner decomposes the adaptive filter into a constrained (ref- 
erence vector filters) and an unconstrained (auxiliary vector filters) 
components. The weight vector is iterated by subtracting the scaling 
auxiliary vector from the reference vector. The scalar factor and the 
auxiliary vector depend on each other and are jointly calculated ac- 
cording to the CCM criterion. The proposed robust AVF algorithm 
provides an iterative exchange of information between the scalar fac- 
tor and the auxiliary vector and thus leads to a fast convergence and 
an improved steady-state performance over the existing techniques. 
Simulations are performed to show the performance and the robust- 
ness of the proposed scheme and algorithm in several scenarios. 

Index Terms — Beamforming, antenna arrays, constrained con- 
stant modulus, auxiliary vector. 

1. INTRODUCTION 

Adaptive beamforming techniques are widely used in numerous ap- 
plications such as radar, wireless communications, and sonar ||I]-||2], 
to detect or improve the reception of a desired signal while suppress- 
ing interference at the output of a sensor array. Currently, the beam- 
formers designed according to the constrained minimum variance 
(CMV) and the constrained constant modulus (CCM) criteria are 
among the most used criteria due to their simplicity and effective- 
ness. The CMV criterion aims to minimize the beamformer output 
power while maintaining the array response on the direction of the 
desired signal. The CCM criterion is a positive measure (Chapter 
6 in 1 2 1) of the deviation of the beamformer output from a constant 
modulus (CM) condition subject to a constraint on the array response 
of the desired signal. Compared with the CMV, the CCM-based 
beamformers exhibit superior performance in many severe scenarios 
(e.g., steering vector mismatch) since the positive measure provides 
more information for parameter estimation with constant modulus 
signals. 

For the design of adaptive beamformers, numerous adaptive fil- 
tering algorithms have been developed using constrained optimiza- 
tion techniques |3|, |5|. The stochastic gradient (SO) and recur- 
sive least squares (RLS) |3|, |5 | are popular methods with different 
tradeoffs between performance and complexity. A major drawback 
is that they require a large number of samples to reach the steady- 
state when the array size is large. In dynamic scenarios, filters with 
many elements usually provide a poor performance in tracking sig- 
nals embedded in interference and noise. The multistage Wiener 
filter (MSWF) 1 8 1 provides a way out of this dilemma. The MSWF 
employs the minimum mean squared error (MMSE) criterion and its 



extended versions with the CMV and the CCM criteria are reported 
in (9), 1101 . Another cost-effective technique is the auxiliary vector 
filtering (AVF) pi 11 algorithm. In this scheme, an auxiliary vector 
is calculated by maximizing the cross correlation between the out- 
puts of the reference vector filter and the previously auxiliary vector 
filters. The weight vector is obtained by subtracting the scaling aux- 
iliary vector from the reference vector. In 11121 . the AVF algorithm 
iteratively generates a sequence of filters that converge to the CMV 
filter with a small number of samples. Its application in adaptive 
beamforming has been reported in [|13l . 

Motivated by the fact that the CCM-based beamformers outper- 
form the CMV ones for the CM signals, we propose an AVF algo- 
rithm based on the CCM design for robust adaptive beamforming. 
The beamformer structure decomposes the adaptive filter into a con- 
strained (reference vector filters) and an unconstrained components 
(auxiliary vector filters). The constrained component is initialized 
with the array response of the desired signal to start the iteration and 
to ensure the CCM constraint, and the auxiliary vector in the uncon- 
strained component can be iterated to meet the CM criterion. The 
weight vector is computed by means of suppressing the scaling un- 
constrained component from the constrained part. The main differ- 
ence from the existing AVF technique is that, in the proposed CCM- 
based algorithm, the auxiliary vector and the scalar factor depend on 
each other and are jointly calculated according to the CM criterion 
(subject to different constraints). The proposed method provides an 
iterative exchange of information between the auxiliary vector and 
the scalar factor and also exploits the information about the CM sig- 
nals, which leads to an improved performance. Simulations exhibit 
the robustness of the proposed method in typical scenarios including 
array mismatches. 

The rest of this paper is organized as follows: we outline a sys- 
tem model and the problem statement in Section 2. The proposed 
scheme is introduced and the CCM-AVF algorithm is developed in 
Section 3. Simulation results are provided and discussed in Section 
4, and conclusions are drawn in Section 5. 

2. SYSTEM MODEL AND CCM BEAMFORMER DESIGN 

Consider q narrowband signals that impinge on a uniform linear ar- 
ray (ULA) of m (m > q) sensor elements. The sources are assumed 
to be in the far field with directions of arrival (DOAs) Oq,. . . ,9q-i. 
The ith received vector x{i) G can be modeled as 

x{i) ^ A{6)s{i) + n{i), i^l,...,N, (1) 

where 6> = [^o, . . . , Og-if G R"""^ is the signal DOAs, A{e) = 
[a{6o), . . . , a{6q-i)] G C'"^'' comprises the signal steering vectors 

^(g^-) ^ [1,6"^"-'^''°'''''=,..., g-2'rj(™-l)^cosefcjT g 

{k = 0, . . . , q — 1), where Ac is the wavelength and d is the inter- 
element distance of the ULA (d = Ac/2 in general), s{i) G C^^ 
is the source data, n(i) G C"*^^ is assumed to be a zero-mean 



spatially white Gaussian process, A'^ is the number of snapshots, and 
(•)^ stands for transpose. To avoid mathematical ambiguities, the 
steering vectors a{9k) are normalized and considered to be linearly 
independents. The output of the beamformer is 

y{i) = w"{i)x{i), (2) 

where — [wi{i), . . . ,Wm.{i)]'^ G C™'^^ is the complex weight 
vector of the beamformer, and (■)^ stands for Hermitian transpose. 

With the signals introduced in ([TJ and ([2}, we can present the 
CCM beamformer design by minimizing the following cost function 

Jcm{w{i)'j = E{ 1^ — !/] subject to {i)a{9o) — 1, 



3.2. Proposed CCM-AVF Algorithm 

From Q, it is necessary to determine the auxiliary vector g^ii) and 
the scalar factor (i) for the calculation of lOi (i). Assuming (i) 
is known, ^i(i) can be obtained by minimizing ¥,{[wi{i)x{i) ~ 
1]^}. Substituting Q into this cost function, computing the gradient 
with respect to (i) and equating it to zero, we obtain 



where is the direction of the signal of interest (SOI) and a{9o) 
denotes the corresponding steering vector. The cost function is the 
expected deviation of the squared modulus of the array output to a 
constant, say — 1. The constraint is set to maintain the power of 
the SOI and to ensure the convexity of the cost function. The weight 
expression obtained from (|3) is 

r^u R-irU (^ [p" ii)R''ii)a{eo) - l]aieo) , 

W{l+1) = R I {p I - i rn\ 1' 

aH{9o)R '-{i)a(9o) 

(4) 

where R{i) =E[\y{i)\^x{i)x"{i)] £ C'"^'",p(j) = E[y* {{)x{i)] € 
C"^^, and (■)* denotes complex conjugate. Note that Q is a func- 
tion of previous values of w{i) (since y{i) = {i)x{i)) and thus 
must be initialized to start the iteration. We keep the time index i in 
R{i) and p{i) for the same reason. The calculation of the weight 
vector is costly due to the matrix inversion. The SG or RLS type 
algorithms can be employed to reduce the computational load but 
suffer from a poor performance when the dimension m is large. 

3. PROPOSED CCM BEAMFORMER DESIGN 
AND AVE ALGORITHM 

In this section, we introduce a CCM-based adaptive filtering struc- 
ture for beamforming and develop an efficient CCM-AVF algorithm 
for robust adaptive beamforming. 



3.1. Proposed CCM Beamformer 

We define the cost function for the beamformer design, which is 



{i))^E{[w"{^x{i)-u]'}, 



(5) 



where x{i) — y* {i)x{i) can be viewed as a new received vector to 
the beamformer and v = lis set in accordance with (|3}. 

To obtain the weight solution for the time index i, we start the 
iteration by initializing the weight vector 

wo{i) = a{9o)/\\a{9o)f. (6) 

Then, we subtract a scaling auxiliary vector (unconstrained com- 
ponent) that is orthogonal to 0(6^0) from (i) (constrained compo- 
nent) and obtain 

lui(i) = lyo(j) - (7) 

where 9^(1) G c""^^ with gf (i)a(6'o) = 0, and fj,i{i) is a scalar 
factor to control the weight of gi{i)- The auxiliary vector is sup- 
posed to capture the signal components mx{i) that are not from the 
direction 6*0. The aim of l|7j is to suppress the disturbance of the 
unconstrained component while maintaining the contribution of the 
SOI. The cost function in (|5} appears in unconstrained form since 
the constraint has been incorporated in the weight adaptation. 



Mi(j) 



g^{i)R{i)wo{{)^g^{i)p{t) 
g?{^)R{i)g,{^) 



(8) 



(3) C 



where R{i) = E[x{i)x" (i)] e C'"""" and = E[x{i)] G 
""^^ Note that = 0, i.e., R{i)wo{i) = p{i) needs to be 



avoided here since the design is equivalent to a matched filter if it 
happens. 

On the other hand, the calculation of the auxiliary vector gi{i) 
should take the conditions gf {i)a{6o) = and g^ (i)g^{i) — 1 
into account. The constrained minimization problem with respect to 
g^ (i) can be transformed by the method of Lagrange multipliers into 



JL{wi{i)) =E{[ii;f(i)i(i)-l]'} 
-2$R{Ai[gf(i)gi(i) 



X2gf{i)a{9o)}, 



(9) 



where Ai and A2 are scalar Lagrange multipliers. For the sake of 
mathematical accuracy, we note that the cost function to be mini- 
mized is phase invariant, namely, if g j (i) satisfies it, so does g^ (i)e-''^ 
for any phase 0. To avoid any ambiguity, we assume that only one 
auxiliary vector can be obtained. 

Following the procedure to get ^1 (i), the auxiliary vector can be 
expressed by 



Mi(«)Py(*) - >^2a(9o) 

a; ' 



(10) 



wherepy(i) = E[(l-y(i))*i(i)] G C'"^^andy(i) = w"{i)x{i). 
We keep the time index i in (i) since it is a function of ■w{i), which 
must be initialized to provide an estimation about y{i) and to start 
the iteration. 

The expression of (i) is utilized to enforce the constraints and 
solve for Ai and A2. Indeed, we have 



Ai 



^lii)a"{9o)p,{i) ^ 

Ml(*)Py(«) u„fO \U9. '^(^O, 



Ao = 



fiUi)a" {9o)pJi) 



(11) 



(12) 



MOoW ' 

where || • || denotes the Euclidean norm. Substitution of Ai and A2 
back in JlOl l leads to g^{i) that satisfies the constraints and mini- 
mizes (with /ii(i)) the squared deviation of y{i) from the CM con- 
dition, yielding 



Ml(*)Py(*) 



lk(9o)ll 



|Mi(i)Py(«) 



Mj(i)a," (So)p,,(*) /o \\ 



(13) 



So far, we have detailed the first iteration of the proposed CCM- 
AVF algorithm for time index i, i.e., wo{i) in (|6}, lyi(i) in 0, /ii (i) 
in ([8}, and gi{i) in ( 113b . respectively. In this procedure, x{i) can be 
viewed as a new received vector that is processed by the adaptive 
filter wi{i) (first estimation of w{i)) to generate the output y{i), in 
which, lyi(i) is determined by minimizing the mean squared error 



Table 1. PROPOSED CCM-AVF ALGORITHM 

For the time index i — 1,2, . . . , N. 

Initialization: 

Iterative procedures: 

Forfc = 1,2 K 



small positive value. 



Mfc-l(i)PyW lk(6)o)||^ °(^0) 



ifSfe(*) = Othen EXIT. 

„ (--s _ gf (i)K(i)«)fc-i(i)-gf (»)p(i) 

Weight expression: 

iy(i) = WK{i)- 



between the output and the desired CM condition. This principle is 
suitable to the following iterations with W2{i),'ws{i), .... 

Now, we consider the iterations one step further and express the 
adaptive filter as 

2 

W2{t) = Wo(i) - ^llk{i)gk{i) = - M2(*)92(*): (14) 

fc=l 

where At2(i) and g2(*) calculated based on the previously 

identified /ii (i) andgj^(i). /i2(i) (m2(*) 7^ O) is chosen to minimize 
the cost function '¥^{[w2 {i)x{i) — 1]^} under the assumption that 
§2 (i) is known beforehand. Thus, we have 



M2(j) = 



i)p{i) 



g^{i)R{i)gS) 



(15) 



The auxiliary vector (*) is calculated by the minimization of the 
cost function subject to the constraints g^(i)a(&o) = and 



9?(*)92(0= 1' which is 



g,ii)= ''^'''^\,.,}^^fLj^'^ . (16) 



Mj(i)a"{eo)Py(») la \ 



The above iterative procedures are taken place at time index i to gen- 
erate a sequence of filters Wk{i) with k = 0, 1, . . . being the itera- 
tion number. Generally, there exists a maximum (or suitable) value 
of k, i.e., fcmax ~ K, that is determined by a certain rule to stop iter- 
ations and achieve satisfactory performance. One simple rule, which 
is adopted in the proposed CCM-AVF algorithm, is to terminate the 
iteration if 9^(1) = is achieved. Alternative and more compli- 
cated selection rules can be found in [13]. Until now, the weight 
solution at time index i can be given by w{i) = WK{i). The pro- 
posed CCM-AVF algorithm for the design of the CCM beamformer 
is summarized in Table [T] 

3.3. Interpretations about Proposed CCM-AVF Algorithm 

There are several points we need to interpret in Table [T] First of all, 
initialization is important to the realization of the proposed method. 

is set to estimate y{i) and so R{i), p{i), and Py{i)- wo{i) is 
for the activation of the weight adaptation. Note that the calculation 
of the scalar factor, e.g., in ((8), is a function of gi{i) and the auxil- 
iary vector obtained from ( 1131 depends on It is necessaiy to 



initialize one of these quantities to start the iteration. We usually set 
a small positive scalar value fJ,o{i) for simplicity. Under this con- 
dition, the subscript of the scalar factor for the calculation of g^, (i) 
should be replaced by fc — 1 instead of k, as shown in Table[T] 

Second, the expected quantities R{i), p{i), and Py{i) are not 
available in practice. We use a sample-average approach to estimate 
them, i.e.. 



I ^ — ' I ^ — ' 

(=1 (=1 



(17) 



where R(i), p{i), and Py(i) are substituted by their estimates in the 
iterative procedure to generate Wk{i)- To improve the estimation 
accuracy, the quantities in ( |17t can be refreshed or further regular- 
ized during the iterations. Specifically, we use Wk{i) in the iteration 
step instead of w{i) in the initialization to generate y{i), and related 

x{i) and y{i), which are employed to update the estimates R{i), 
p{i), and Py{i)- Compared with w{i) = a{6o)/\\a{6o)\\^ , Wk{i) is 
more efficient to evaluate the desired signal. Thus, the refreshment 
of the estimates based on the current Wk{i) is valuable to calculate 
the subsequent scalar factor and the auxiliary vector. 

Third, we drop the normalization of the auxiliary vector (121- 
1 14|. Note that the calculated auxiliary vectors g^ii) are constrained 
to be orthogonal to a {60). The orthogonality among the auxiliary 
vectors is not imposed. Actually, the successive auxiliary vectors 
do satisfy the orthogonality as verifies in our numerical results. We 
omit the analysis about this characteristic here considering the paper 
length. 

The proposed CCM-AVF beamformer efficiently measures the 
expected deviation of the beamformer output from the CM condi- 
tion and provide useful information for the proposed algorithm for 
dealing with parameter estimation in many severe scenarios includ- 
ing low signal-to-noise ratio (SNR) or steering vector mismatch. The 
proposed CCM-AVF algorithm employs an iterative procedure to ad- 
just the weight vector for each time instant. The matrix inversion 
appeared in (|4j is avoided and thus the computational cost is lim- 
ited. Since the scalar factor and the auxiliary vector depend on each 
other, the proposed algorithm provides an iterative exchange of in- 
formation between them, which are jointly employed to update the 
weight vector. This scheme leads to an improved convergence and 
the steady-state performance that will be shown in the simulations. 



4. SIMULATIONS 

Simulations are performed for a ULA containing m = 40 sensor el- 
ements with half-wavelength interelement spacing. We compare the 
proposed algorithm (CCM-AVF) with the SG 13], RLS |5|, MSWF 
|_10|, and AVF 1 12| methods. With respect to each method, we con- 
sider the CMV and the CCM criteria for beamforming. A total of 
1000 runs are used to get the curves. In all experiments, BPSK 
sources' powers (desired user and interferers) are cr| — af = 1 and 
the input SNR= dB with spatially and temporally white Gaussian 
noise. 

Fig. [T]includes two experiments. There are g = 10 users, in- 
cluding one desired user in the system. The scalar factor is /io (i) — 
0.01 and the iteration number is K = 3. In Fig[T] (a), the exact 
DOA of the SOI is known at the receiver. All output SINR values 



increase to the steady-state as the increase of the snapshots (time in- 
dex). The RLS-type algorithms enjoy faster convergence and better 
steady-state performance than the SG-type methods. The proposed 
CCM-AVF algorithm converges rapidly and reaches the steady-state 
with superior performance. The CCM-based MSWF technique with 
the RLS implementation has comparative fast convergence rate but 
the steady-state performance is not better than the proposed method. 
In Fig. [T](b), we set the DOA of the SOI estimated by the receiver to 
be 1° away from the actual direction. It indicates that the mismatch 
induces performance degradation to all the analyzed algorithms. The 
CCM-based methods are more robust to this scenario than the CMV- 
based ones. The proposed CCM-AVF algorithm has faster conver- 
gence and better steady-state performance than the other analyzed 
methods. 




(a) Number of Snapshots (b) Number of Snapshots 

Fig. 1. Output SINR versus the number of snapshots for (a) 
ideal steering vector; (b) steering vector mismatch 1°. 

In Fig. |2l we keep the same scenario as that in Fig. [T](a) and 
check the iteration number for the existing and proposed methods. 
The number of snapshots is fixed to TV = 500. The most adequate 
iteration number for the proposed CCM-AVF algorithm is if = 3, 
which is comparatively lower than other analyzed algorithms, but 
reach the preferable performance. We also checked that this value is 
rather insensitive to the number of users in the system, to the number 
of sensor elements, and work efficiently for the studied scenarios. 



5. CONCLUDING REMARKS 

We developed an AVF algorithm based on the CCM design for robust 
adaptive beamforming. The algorithm provides a positive measure 
of the expected deviation of the beamformer output from the CM 
condition and thus is robust against the severe scenarios. The weight 
solution is iterated by jointly calculating the auxiliary vector and the 
scalar factor, which iteratively exchange information between each 
other and lead to an improved performance over prior art. The selec- 
tion of the iteration number may be more efficient and adaptive with 
the change of the system (e.g., the number of users change) if we em- 
ploy other techniques |13|. We will consider further improvements 
to the proposed CCM-AVF algorithm in the near future. 
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Fig. 2. Output SINR versus the number of iterations. 
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